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Abstract 

We work in a model where all CP violating phenomena have a common source. CP is 
spontaneously broken at a large scale V through the phase of a complex singlet scalar. An 
additional SU(2)l singlet vector-like down-type quark relates this high scale CP violation 
to low energy. We quantitatively analyze this model in the quark sector. We obtain the 
numerical values of the parameters of the Lagrangian in the quark sector for a specific 
ansatz of the 4x4 down-type quark mass matrix where the weak phase is generated 
minimally. Zbb vertex will modify in presence of the extra vector-like down-type quark. 
From the experimental lower bound of the partial decay width Z — > bb we find out the 
lower bound of the additional down-type quark mass. Tree level flavor changing neutral 
current appears in this model due to the presence of the extra vector-like down-type quark. 
We give the range of values of the mass splitting Am B q in B® — B® system using SM box, 
Z mediating tree level and Z mediating one loop diagrams together for both q = d,s. 
We find out the analytical expression for T\ 2 in this model from standard box, Z and 
Higgs mediated penguin diagrams for — B® system, q = d,s. From this we numerically 
evaluate the decay width difference | Ar s q /^B q \ ■ We also find out the numerical values of 
the CP asymmetry parameters aj and for the decays B^ — > J/tpK s and B® — > 7r + 7r~ 
respectively. We get the lower bound of the scale V through the upper bound of the 
strong CP phase. 



1 Introduction 



CP violation is an important phenomena in the context of particle physics and cosmology. CP 
violation is directly observed in the decays of K and B mesons. The present experimental re- 
sults d] are consistent with the standard model(SM). The single phase present in the Cabibbo, 
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Kobayashi, Maskawa(CKM) mixing matrix is responsible for this kind of CP violating phe- 
nomena. Present observational data of neutrino oscillations prove that neutrinos are massive. 
Minimal extension of the SM with neutrino mass generate neutrino mixing similar to the quark 
mixing. So the phases in the neutrino mixing matrix generate CP violating phenomena in 
the leptonic sector. The strong CP violation comes from non-perturbative instanton effects in 
the SM. This leads to the so called strong CP problem for which various solutions have been 
proposed. The bound on the electric dipole moment of neutron gives the bound on strong CP 
phase. In the context of cosmology baryon asymmetry in the Universe (BAU) gives an obser- 
vational evidence for CP violation. Decay of heavy Majorana neutrino to lepton(both charged 
and neutral) and scalar(both charged and neutral) generate lepton asymmetry which violates 
CP P]. One way of generation of baryon asymmetry is through sphaleron mediated process 
from lepton asymmetry. 

Different CP violating processes are not related in general. So the question may arise 
whether it is possible to find a model where all kinds of CP violation have a common origin. 
In fact there are a few models [31 H] . We work in the model proposed by Branco, Parada and 
Rebelo(BPR) |3j. CP breaks in this model spontaneously through a SU(3) C x SU(2) L x U{l) Y 
singlet complex scalar which gets a vacuum expectation value with nonzero phase, ^exp(zct). 
This large scale phase is responsible for all kinds of CP violation. A vector-like down-type quark 
and three generations of right handed neutrinos relate the large scale phase to low energy CP 
violating phase in the quark sector and the leptonic sector respectively. We quantitatively 
analyze this model for the quark sector. In Section |2] we describe the model [3] in detail. 
In Section |3] we calculate parameters of the Lagrangian for a specific ansatz for the down- 
type quark mass matrix, under the assumption that additional down-type quark mass is larger 
than the standard down-type quark masses and using the experimental range of values of the 
CKM parameters and the standard down-type quark masses as the inputs. We take CKM 
matrix in "standard" parametrization j3] for convenience. Presence of a vector-like down-type 
quark generate flavor changing neutral current (FCNC) at the tree level. It changes the flavor 
preserving vertex as well. 

We see the effect in Z — > bb decay in Section El In this Section we get the lower bound 
of the additional down-type quark mass from the experimental lower bound of Z — > bb decay 
width. We also see the effect in Amg for q = d, s of this additional down-type quark mass in 
the same Section El We also discuss the effect of this model on Ar# in the same section for 
both q — d,s. We evaluate the expression of linearly Um dependent contribution to T\ 2 which 
is require to calculate the decay width differences. Um is tree level coupling of b and d quarks 
with Z and also with neutral physical as well as fictitious scalar. We only consider the light 
Higgs contribution and disregard heavy scalar contribution. At last in this section we calculate 
the numerical values of the CP asymmetry parameters aj and in the decays B% — > J/ipK s 
and B% — > tt + tt~ respectively. Using the parametric solution and the extra down-type quark 
mass bound to strong CP phase we get the lower bound on the scale V in Section |U Section 
contains concluding remarks. 
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2 The Model 



The particle content of the standard model is 

9l=(^)(3,2^), <4(3,1^), 4(3,1,-1), 

/l=(^J(1,2,~), 4(1,1,-1), *=( ) (1,2,|). (1) 

q' L , ip l ' L and $ are the standard quark, lepton and scalar doublets respectively. Here u' L , 
u' R , d' L and d' R are standard three generation left handed up-type, right handed up-type, left 
handed down- type and right handed down- type quark fields respectively l' L , l' R and v' L are 
standard three generation left handed lepton, right handed lepton and left handed neutrino 
fields respectively 

The model of BPR [3] has the same gauge symmetry as the SM, viz., SU(3)c x SU{2) L x 
U(1)y- But in addition it has discrete symmetries Z 4 x CP. There are also six additional 
multiplets in this model along with the SM particles. One is a totally gauge singlet complex 
scalar S. Others are the two chiralities of an SU(2)l singlet vector-like down-type quark D' , 
and three generation of gauge singlet right handed neutrinos N' R . These additional multiplets 
have the following SU(3)c x SU{2) L x U(l)y gauge group representations 

^(3,1,-1), 2^(3,1, ~), A^(1,1,0), 5(1,1,0). (2) 

The primes on the fermion fields in Eq. and Eq. (j2J) imply that they do not necessarily have 
definite mass. Under Z± the fields transform as, 

S^-S, D'^-D', v'i-n^, Ir-+U'r, K^iK (3) 

while all the other fields remain invariant. The Yukawa interactions in this model are, 

- Jgf = q' L Y d $d! R + q' L Y^u' R + (f q S + f q 'S*)D' L d' R + ^D' L D' R + h.c. (4) 

in the quark sector and 

= fiY l ^l' R + ^M' R + N' T R C~\US + fJS*)N' R + h.c. (5) 

in the leptonic sector, where $ = ia^* and S* is the conjugate field of S. Y u , Y d in Eq. (J3J) 
and Y l in Eq. (jSJ) are the standard 3x3 Yukawa coupling matrices. f q , f in Eq. (jjj) are the 
nonstandard 1x3 Yukawa coupling matrices. Y u , f u and f' v in Eq. (jSj) are also nonstandard 
3x3 Yukawa coupling matrices. Without loss of generality we choose the basis of the fields 
such that Y u is diagonal with the real positive entities. C in Eq. (0) is a 4 x 4 matrix which 
relate any charge conjugated four spinor field to its original field, ip c = 'joCip*. The Lagrangian 
Jzf ? y in Eq. (jlj has bare mass term of D' where [Iq is the bare mass parameter. Presence of this 
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term does not violate any symmetry of the Lagrangian. All the couplings of the Lagrangian of 
this model are real due to CP symmetry. The most general gauge, CP and Z 4 invariant scalar 
potential [6 of this model is 

y = _ At 2 $t$ + A($t$) 2 + SS^^ + bxSS*) + {S 2 + S* 2 ){a 2 + b 2 SS*) 

+h{S* + S* 4 ) + ($ t $){c 1 (5 2 + S* 2 ) + c 2 SS*}. (6) 

Because of the terms a 2 ,b 2 , &3 and c\ the field S can acquire a complex vacuum expectation value 
(VEV) : (S) = ^exp(ia). This VEV breaks CP and Z A spontaneously 1 . The SU(2) L x U{l) Y 

symmetry spontaneously breaks when the scalar $ takes VEV, (<&) = ( „ V It is natural to 
think that V ^> v due to the fact that the scalar S is a gauge singlet. 



3 Quantitative analysis in the hadronic sector 



The VEV of $ and S break the SU(3) C x SU(2) L x U(1) Y x Z 4 x CP symmetry down to 
SU(3)c x U(1)q x Z 2 . This generates mass terms. Mass part of the Lagrangian in the quark 
sector is, 



JS7 = d' L m° d d' R + u' L m° u u' R + M° D D' L d' R + ^D' L D' R + h.c. (7) 



where, 



and 



M° D = ^={f q ex V (ia) + //exp(-ia)}. (9) 

Since we are in a basis where Y u is diagonal with the real positive entities, the up-type quark 
mass matrix m" will also be diagonal, m" = diag(m u , m c , m t ). So the up-type fields v! are 
physical. We will call them u(unprimed) for future reference. This is our convention. The 
down-type quark mass term from Eq. (JJJ) can be written by the following form 

V \ ( W>3X3 {0} 3 xl \ f d' R x 



-sr = l4 ^){{^ } Z wrjl^J' (10) 

The above mentioned 4x4 down-type quark mass matrix is denoted by ^ . jtft can be 
diagonalized by bi-unitary transformation, 



1 Spontaneous breaking of any discrete symmetry generates cosmological domain wall problem. Solution of 
this problem goes beyond the scope of this article. 
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where = diagfrn^, m s , m.&) and Md is the mass of additional down- type quark. Let us 
consider the forms of and as follows 

* = ( K q 3xs ^ 3xl ), *'=( K q h 5 xi ) ■ (12) 

\ ^1x3 ilxl / \ ^1x3 J lxl / 

Relations of the physical basis(unprimed) of the down-type quark fields to the original basis 
are 



d L \ ( Kl & \( d' L 
D L \ B) T* { D'i 



and 



( d R \ = ( K« S' ] \ ( d' R 
\D R ) { T'* J { D' R 

Clearly % diagonalizes 



(13) 



(14) 



M 2 D 



(15) 



This gives the approximate relation [3J E] 

^ = Krh 2 d K- 1 (16) 

where 

^^x'- '-^y' (it) 

is the effective mass matrix square for the standard down-type quarks under the assumption 

f « 1. («> 

and with 

Ml « M°M° f + ,4. (19) 

X is the CKM matrix and it is approximately unitary for the hermiticity of M 3 . The other 
blocks of % in Eq. (fT2|). under the assumption in Eq. (JTSJ) obtain the forms [S| 

fiK ibf- Tssl (20) 

and 



The constraints of unitarity of the matrix % defined in Eq. (I12|) . imply the relations 

Y,\Ki^ = l-\Ri\ 2 (22) 

3 

for all generations i, i — u, c, t. This implies that the sum of all couplings of any of the up-type 
quarks to all the standard down type quarks is not same for all generations. Presence of \Ri\ 
term in the right hand side of Eq. (|22jl breaks the weak universality. Smallness of the 2nd term 
in the right hand side of Eq. (f2*2*j) due to the assumption in Eq. (fTHj) suggests that unitarity of 
CKM matrix K and weak universality hold only approximately in this model. 



Here the CKM phase can be generated through the second term of M' in Eq. ([17)1 . Note 
that this term is not suppressed by the scale V because and Mp, as we have seen from Eq. 
(jHJ) and Eq. (fTTIJ) . both have dependence on the scale V. 

Now we try to analyze this model quantitatively. Since up-type quark mass matrix is 
diagonal, the remaining parts of =$f g m in Eq. (J7J have 18 real parameters. These parameters are 
9 of m", 3 of f q , 3 of /' and other 3 are a, V and /xo- The matrix M in Eq. (fTTj) and hence 
3^ in Eq. (fTTj) contains all these parameters. The question may arise whether it is possible to 
find realistic solution of the Eq. (jlfijl for the entire range of values of the CKM parameters and 
the masses of the standard down-type quarks. Since ffl is hermitian, there are 9 independent 
coupled equations consisting of these 18 parameters as variables. The reduction of the number 
of variables can be obtained as the following way. First we assume that m d is symmetric. Then 
we note that and occur in Eq. (fTTj) in the combination 



M, 



v 



(23) 



which can be parametrized as follows 



F 1 



( /iexp(za) + exp( 

f 2 exp(ia) + ftexp(-ia) |, (24) 
V /sexp(ia) + ftexpi 



-ia 



where fi — :(f q )i and f- = 

y/5M D (fq)i- Now if we consider f — f, the matrix Jif in Eq. ((TTJ 
becomes real symmetric. This implies that the CKM phase is zero. Now let us see minimally 
how we can generate a non zero CKM phase. To this end, we take the ansatz 

h = f[ = f 
h = & = Pf 

/s = / + A/ 

fs = /■ (25) 

So ultimately we have 10 variables which are /, Af, a, (3 and six elements of m? d . So we are in 
a position where we can find the parameters of the Lagrangian from Eq. (|16|). We choose the 
values of the CKM parameters and the masses of the standard down-type quarks from their 
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1 T r~Y~~\ Q r\T mT\llf 1— ' o V O "IT! O 1 QVC 

lypc Ol HipUL JT dldlilcLcIo 


opeciiicanon 


n/xpei lmeniai rvange 


rtaiiciom v^noice 


CKM Parameters 


sin ^12 


n 9997 fn n 99^Q 


n 99^8 


sin ^13 


u.uuoz to u.uu^tz 


n nn f ?7 


sin ^23 


n d^qs tn n n498 

U.UOc/O LO U.U^ZO 


n n^QS 
u.uoyo 




46° to 74° 


73.6818° 


Mass Parameters 


m d 


4.0 to 8.0 MeV 


5.9171 MeV 


m s 


80.0 to 130.0 MeV 


81.6310 Mel/ 


m b 


4.1 to 4.4 GeV 


4.3358 GeV 



Table 1: Here fly's are the CKM angles and 5c km is CKM phase in "standard" parametrization 
of CKM matrix [5j. Experimental range of the sin of the CKM angles, CKM phase and the 
standard down-type quark masses are taken from PDG [I]. In [T] the values quoted for mass 
of md, m s are measured at the scale 2 GeV and = fh^fhb) in the MS scheme. We randomly 
choose a set of the above inputs from their experimental range. 



experimental range as in Table d . Using these values in the right hand side of the Eq. (jl6|) 
we solve this equation for different values of (3. We observe that there is a range of the values 
of f3, —3.2 < (3 < 2.12, where Eq. (jl6j) fails to give real solutions for the parameters of ffl for 
the entire range of values of the CKM parameters and the standard down-type quarks masses 
as in Table [TJ We get the real solutions outside this range of values of (3. 

To get the essence of the solutions let us show that how we have proceeded to reach the 
goal. First we write m d and Fjj in the following forms: 

/ Xi x 2 x 3 \ ( 1 \ 

m d = v x 2 x A x 5 F D T = ^fx>j (3 . (26) 

\ x% x b x & ) \ x 8 + ix 9 J 

Here x\ to x§ are six independent dimensionless mass parameters of m° d . x-j to xg are three 
independent dimensionless parameters of Fp. We have written Fp in Eq. (|26j) with a pre-f actor 
y/xj. This makes x-j dependence of Eq. (Ilfij) linear. The xj to xg are related to /, Af and a in 
the following way 

,,22 , . A / A/ tan a 
x 7 = 4/ cos a, x 8 = 1 + — , xg = — . (27) 

For every set of inputs we have four sets of solutions. These can be divided into two pairs of 
sets, where within a pair, the two solutions are related only by change of signs for x\ to xq. 
The solutions of different pairs have the values of the same order of magnitudes. We provide 
two distinct sets of solutions in Table 121 for a set of inputs in the rightmost column of Table 
Qwith (3 = 3.5 and (3 = —3.5. Inputs are chosen randomly from their experimental range in 
Tabled 

Now let us see that what are the outcomes of these solutions. m d becomes explicitly known. 
Fd also becomes known meaning that Yukawa parameters related to the coupling of the stan- 
dard down-type quarks, extra vector like down-type quark and the singlet scalar also become 
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V ClI IdUlCo 


Results For /3=3.5 


Results For (3=— 3.5 


OCLj INC. -L 


kJCL INC. 




Qpf M n 

kJCL 1NU. £ 


T-i 


n nnm 6?84i 


0001 68^6 


—0 000008467 


—0 00001 ^441 3 


Tn 


n nnrnm 087 


000S90S8Q 


0009Q84Q1 


000^?9Q01 


X3 




00099Q067 


— 000^904QQ 


0009Q1 798 




—0 000122682 


—0 0000551233 


000131765 


0000463298 




-0.000845869 


0.000848642 


0.00107484 


-0.00106796 




-0.0212162 


0.0212857 


0.026958 


-0.0267854 




0.0488582 


0.0539958 


0.0289893 


0.0353266 


Xg 


-2.53447 


2.14896 


4.41581 


-3.72177 


x 9 


-0.849496 


0.766168 


-1.28386 


1.06035 



Table 2: Two distinct sets of solutions of the Eq. (|T6|) for the inputs of the rightmost column 
of Tabled with P = ±3.5. 



Quantities 


Results For /?=3.5 


Results For /3=— 3.5 


Set No. 1 


Set No. 2 


Set No. 1 


Set No. 2 


/ 


±0.113667 


±0.139648 


±0.0909458 


±0.0963174 


A/ 


T0.803504 


±0.3209 


±0.621307 


^0.909577 


a 


13.5145 


33.6967 


-20.5991 


-12.6567 




0.996472 


0.996495 


0.997163 


0.997127 



Table 3: Values of /, A/, a and F D F* D for the solution in Table |21 using Eq. (J57J). 



known upto a factor Mjj/V. Using the relations in Eq. (}2"Tj) we obtain the values of /, Af 
and a in Table El for the different solutions in Table 121 The R block of as in Eq. ()12j) is 
the coupling of charged current interactions among up- type quarks u and D. S^S is a matrix 
which generates FCNC for the standard down-type quarks. Using Eq. (J23|) . we can rewrite Eq. 
(j2Ujl and Eq. (J2U as 
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M. 



D 



M. 



(28) 



The numerator of R and S in above Eq. ()28|) become known from the numerical solution. Their 
forms for the Set-1 solution with f3 = 3.5 in Tableware 



R 



i 

Mn 



/ 0.101309 -^0.0117825 
0.109596 -^0.0390722 
V 2.74902 -^0.980013 



S T 



i 



-0.0993492 - z0.000938182 
-0.0228038 - iO.000215232 
-2.75133 - z0.980436 



(29) 



Now let us observe two parts of the M 3 matrix in Eq. (|17|). Using the Set-1 solution for (3 = 3.5 
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in Table |21 we have 



m° d m d 



0.0110282 
0.0137888 
0.309574 



0.0137888 
0.0496956 
1.08766 



0.309574 
1.08766 
27.287 



(30) 



and 



m° d M D ] M D m° d ] 

M " 
0.0104023 
0.0115634 -i0.00266705 
0.290048 - i0.0668936 



0.0115634 + i0.00266705 
0.0135379 
0.339572 + ^5.14869 x 10~ 6 



0.290048 + i0.0668936 
0.339572 -z5. 14869 x 10" 6 
8.51754 



(31) 



It is interesting to see that two parts of ffl contributes roughly equally for many elements. 
The second part of ^ in Eq. (}3*Tj) is new physics term and it is not suppressed. This feature 
is independent of the inputs we are giving. Another input independent feature is that the 
contribution of bare part to the mass of D quark Mb in Eq. ()19j) is smaller than the Yukawa 
part. The Yukawa part M° D M D ] / ' M 2 D = F D F D contributes approximately 99.64% to M D as in 
Table El for the solutions in Table 2 with (3 = 3.5. The contribution for the solutions (3 = —3.5 
is nearly 99.71% as in Table El Although bare mass term is gauge invariant, its contribution to 
the mass of D quark is subdominant. One conclusion we may draw from here that Mb should 
be less than V for the validity of perturbation theory. Now let us see that how we can get more 
information about Mb from Z — > bb decay and how Mb value affects — B d mixing. 

The presence of extra vector quark generates FCNC at the tree level. It also changes the 
flavor preserving neutral current. The interaction of Z and the standard down- type quark is 



C/>d 



where d\ = d, d 2 = s, d 3 = b, 



and 



where 
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2 cos 9 W £-± x 



(9Ah = gTkj - m 



■1 j 



and the tree level couplings in SM are 



9 S v M 
9 S A M 



u = s ] s 



J 3 d -2Qsin 2 6 W , 



(32) 



(33) 



(34) 



(35) 



(36) 



9 



with J| = —1/2, Q = —1/3. The flavor preserving vertex appreciably changes by loop cor- 
rections. We want to see the effect of loop correction along with the new contribution of U bb 
dependent terms in Z —>■ bb process [7j. The partial decay width of Z — > bb in terms of effective 
coupling constants g Vcff and g b A is [8J 

= [(i - ±y)k\g b v eS l 2 (i + 2y) + \g b A J 2 (l - %)} x (l + 5 QED ) 

H\9vJ + \9Aj 2 )kc D ] + A QCD (y) (37) 
where y = ra b /M z , QED correction Sqed and factorizable QCD correction Sqcd are 
$qed = 3aQ 2 /4:7i, 

Sqcd = (a. fir) + lAl(a s /n) 2 - 12.8(a s /n) 3 - aa s Q 2 /47r\ (38) 

and the effective couplings are 

9v cS = ^p- b (I b 3 ~Ip bb -2Qsm 2 6 W K b ) (39) 
9a cH = sfp h {ll-llU bh ). 

Using general expression for p b and K b [HI EJ OH] with three loop QCD correction and two loop 
electroweak correction, we obtain 0.9935 < p b < 0.9941 and 1.0341 < n b < 1.0382 for SM using 
experimental range of values of a s (Mz), m t , Mw, Mz and m b p. Hence, we get that the SM 
prediction is 0.37655 < V z ^ h i < 0.37869 GeV. Aqcd(u) in Eq- (|3*7j) is b quark mass dependent 
QCD correction. We disregard here small non-factorizable QCD correction. New physics in 
the form of U bb can decrease g\r , g\ and hence the decay width T z ^ b i, since by definition 
< U bb < 1. Experimental dataat the la level given 0.37593 < T z ^ b ~ b < 0.37912 GeV pQ. If 
we take the SM contribution at the theoretical lower limit 0.37655 GeV, we would then need 

U bb < 7.204 x 10" 4 (40) 

in order that the new physics contributions do not violate the experimental bound. From 
the definition of matrix U in Eq. (|3*3j) and using the form of S in Eq. f!28|) we see that U is 
proportional to 1/M|, which is unknown. Hence we can use the upper bound on U bb to obtain 
a lower bound on Mr,. We provide these bounds in TableEJfor different values of f3. The upper 
bound of the value of the other elements of U become known from the lower bound on Mo- 
For the different sets of solution Mp values are approximately same for a given set of inputs. 
The least value of D quark mass Mb ~ 40 GeV for the given range of (3 in the Table 0] The 
maximum value among the elements in matrix is b quark mass which is nearly 4 GeV. So 
the assumption m^/M|,<C 1 in Eq. (|18|) is not too bad at all. This assumption is better for 
smaller f3 region where Mp has to be larger. 

We now want to see the effect in flavor changing process. U dependent extra piece in gy 
and g A in Eq. ()33|) and Eq. ()34j) are responsible for the FCNC. These generate tree level Z 
mediating B® — B® mixing, where q = d,s. The experimental value of the mass splitting Am^ 
in B G d - B° d system is well known, Amf = (3.304 ± .046) x 10~ 10 MeV p. For B° s - B° s system 
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only lower bound exists, Am Bs > 94.8 x 10" 10 MeV at 95% CL p. The contribution of SM box 
[llj, Z mediated tree diagrams ^21 E3] and Z mediated one loop diagrams together in Arris 
has been explictly calculated by Barenboim and Botella JI] 2 , who give off diagonal term of 
B® — B® mixing matrix M\ 2 as 



Ml 



12 



Mll M A 



bq 



where 



M- 



qSM 
12 



G 2 F M^ VB m Bq (B Bq p Bq 
12vr 2 



and 



q*\2 



(41) 



(42) 



Here a and b are 



A bg =l- a(U bq /a) - Ku b q/a 



, C(x t ) 



where 



E(x t ) 
C{xt) 



q\2 



(43) 



E{x t ) 
in sin 2 9 W 



a 



E(x t ) 



(44) 



— Ax t + Hx 2 — x 



4(1 -a:*) 1 



3xj In x t 



2 (l-x t ) 



4 



1 - x t 



x t 3x t \nx t 



x t, 



(45) 



The parameters involved in M\ 2 M in Eq. (|42|) and f Bq are mass and decay constant of 
B® meson respectively. B B is Renormalization Group invariant parameter. i] B is QCD factor. 



Its value is nearly 0.55 [^j. fn Bd = 5279.4 ± 0.5 MeV and m Bg 



x t = rrit/Myy £f = V t lV t q, q — d,s. The uncertainty in the calculation come from 



5369.6 ± 2.4 MeV 0. 

BbJb„- 



B B J Bd = 221 ± 28_t 2 l 2 MeV and s/B~b~Jb s = 255 ± 31 MeV which are obtained from the 
lattice QCD calculations in ^B]. Mass splitting Am Bq is defined in terms of M\ 2 



Am Ba 



2\M? 2 \ 



(46) 



The range of SM predictions for the range of values of m 4 , Mw, triBq and ^jB Bq f Bq are Am|J f 
(2.292 to 5.318) x lO" 10 MeV and Amff = (81.991 to 146.373) x 10~ 10 MeV. Here CKM 
matrix element V t d, V ts and V t b are same as K$i, K32 and K33 respectively. We use V t d = 



0.0079 - i0.00347, V ts 



-0.03906 - i0.000796 and V tb = 0.9992 which are obtained using 



randomly chosen CKM parameters in the rightmost column in Table [0 The range of Am Bs 

2 Although the authors in calculated Amg for B® — system, it can be easily extended to B® — B a s 
system also. 
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and Am^ for the different values of (3 are shown in Table 0] and respectively. We see that 
the range has a small negative shift in the negative (3 region and a small positive shift in the 
positive (3 region from SM range for B® — B® system. The feature is reversed for B® — B° s system 
where range of ArriB a has a small positive shift in the negative (3 region and a small negative 
shift in the positive (3 region from SM range. 

Now we want to see the new physics effect on the decay width difference Ar# 9 in B® — B® 
system for both q — d,s. The SM results come only from box diagrams Presence of tree 
level FCNC coupled to Z and physical Higgs generate new contributions to the decay width 
difference in the B® — B® system through penguin diagrams. We disregard Higgs mediated 
contributions to Mf 2 which are subdominant compared to the Z mediated contributions. We 
now want to calculate the absorptive part of the amplitude r* 2 for B® — > B® in this model, 
keeping only terms linear in Ub q . Standard box, Z and Higgs mediated one loop penguin 
diagrams give 



1 12 



G 2 F m 2 b m Bq (B Bq fl 



2 ■» 

q-> B q / 



8tt 



f$M + u bq J2tf {2(/box(*, i) - /box(i, i)) + /, 



pen \ 



/pen(l)} 



i=2 



(47) 



where SM result due to box diagrams come through fs M (z c ) [HJIH] whose expression is 



fhAzc) = a 2 + \mi{z c + \z\ - \zi) + ef {vr^a 



-z c ) + -z c + -z. 



-zi - r 



(48) 



The absorptive function ^T] related to SM box is 

1 



(1 + ^)(3 - (zi + Zj ) - ( Zi - Zj ) 2 ) \ + 2x b ( Zi + Zj ){zi + Zj - 1) 



(49) 



whereas the absorptive function related to penguin diagrams mediated by Z and light Higgs is 



/pen(i) = 4 x 



3(1 - Xi)(x b - x z ) 



Zi{9l - 9l){l - J - f ) + f (1 + 2* - \xi) + ^(1 + 2z t 



-| \ 9 4 (1 + 2zi) + ^(1 - ^ + ^-(1 - x t ) + + 2z l ){ X -^^) + 

8 13 3 Ixz oxz x b — xh ^xz x b — xh 



Here & = V* b V iq , z, = m 2 /m 2 , x t = mj/M 2 ,, x b = ml/M*,, x z = M'/M 2 ,, g u L = 1/2 - 
2/3sm 2 9w and = — 2/3 sin 2 6 W . i,j = (1,2,3 = u,c,t). The expression /box(^j) an d fpen{i) 
are only for i,j = 1,2. Presence of top inside the loop does not contribute to the absorptive 
part due to the kinematical impossibility. For penguin diagrams with single up-type quarks 
inside the loop does not contribute to the absorptive part. Here the function in Eq. (fKUJ) is for 
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two up- type quarks inside the loop. The 4 factor present infront of the expression of f P en{i) m 
Eq. (JSUj) is due to the fact that there are four types of diagrams contributing equally. These 
are S and T channel diagrams with loop at the different vertices. Now the definition of AT B 

Id is 



B a 



4Re(M 12 rf 2 ) 



Ar 

We numerically evaluate the range of 
constants and life times of B®. m c = 

(1.536 ±0.014) x 10~ 12 s and r Bs = l/T Bs = (1.461 ± 0.057) x 10" 12 s 0. 
the other parameters are given previously. The SM range of AT Bd /T B 

and 



(51) 



Ar B vr Bi 



AT Bq /r B q using range of quark masses, B® masses, decay 
1.15 to 1.35 GeVpQ where m c = m c {m c ). r Bd = l/r_B d = 

Experimental value of 
= (0.3845 to 0.9066)% 

AT B "/T Bd < 18% 



v 12.5191 to 24.0086)%. The exprimental constraints are 

< 54% at 95% CL [HUH]- In our analysis we disregard the 
masses of d, s and u quarks. We also disregard the contributions from heavy scalar meidiated 
penguins to T\ 2 in Eq. (14 7|) . Their contributions are suppressed by their mass. We have shown 



at 90% CL [HEU and AT B '/T B . 



our results of 
parameter (3. 



AT B <'/T Bd in Table Eland 



We have seen that there is 



in Table |U for the different values of the 



AT B °/T Bs 

small positive shift of decay width in the negative 
(3 region and small negative shift in the positive f3 region for the B® — B^ system. On the 
contrary for B° s — B® system we have seen that there is small negative shift of decay width in 
the negative (3 region and small positive shift in the positive (3 region. 

Single vector like down type quark model has simple expression for CP asymmetry of B® 
decay to two channels J/ipK s and 7r + 7r~. Those are 



aj 



T(B° -> J/j)K a ) — T(B° — » J/j>K a ) 
sm{Am Bd t){T(B° -> Jj^K s ) + T(B° -> J/i/>K a )) 
sin(2/3 - arg(A bd )) 



(52) 



and 



T{B G 



T(B° 



IX IX 



sm(Am Bd t)(T(B° -> tx+tx-) + T(B° -> tx+tx-)) 
sin(2a + arg(A M )) 



where a and j3 are defined as usual 



a = arg(--|) 



(53) 



(54) 



(3 = arg(-||). (55) 

We have calculated the numerical values of those CP asymmetry parameters aj and a n in 
this model. The SM value of those parameters are aj = 0.7370 and a n = 0.2549 using CKM 
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elements with randomly chosen CKM parameters in Table ^ The quantity is defined in 
Eq. (|4~5j). We disregard the variaton of arg(A M ) with respect to the change of parameters inside 
for a particular (3. This variation is smaller compared to the variation with respect to (3. 
So, inspite of getting range we have a value for those CP asymmetry parameters for a particular 
(3 value. We show the results in Table El We have seen that value of the asymmetry parameters 
have been changed in the third place and somewhere in the second place after decimal for 
different values of (3 from the SM values. 

Here we cannot have any information about the scale V because of the cancellations of 
the scale from the numerator and the denominator in the second part of Jif as Yukawa part 
dominates over bare part in the mass of D. So the solutions are scale independent. To know 
about the scale we look for use of these solutions in strong CP. 



4 Constraints from the Strong CP Phase 

The existence of topologically nontrivial gauge transformations, and of field configurations 
which make transitions between the different topological sectors of the theory, leads to the 
existence of the new term Q-^^F^ u F a ^ 11 in QCD action. This term violates parity P and 
time-reversal T. It will violate CP due to CPT invariance. Diagonalization of the mass 
matrix requires different U(l) rotation of the left handed fields and the right handed fields. 
So it effectively generates nonzero chiral rotation. This modifies the 9 to 6, where = 
+ arg(detm), where m is quark mass matrix which appear in ^/lui^r, where *Sf contains 
all quarks. Now the model in which we are working, CP is invariant in the Lagrangian and 
breaks only spontaneously. So we cannot keep the term in our Lagrangian. The mass terms 
are generated through spontaneous breaking of SU{2) L x U(l)y and Z 4 symmetry. Since is 
zero, we can write = arg(detm). At the tree level det m is real and hence is zero. We 
should look for one loop correction of the quark mass matrix. Now, loop correction to is 

50 = arg{det (m - £)} 

= -ImfTrtin^E)}, (56) 

where S is self energy matrix. The one loop correction to the mass matrix has been explicitly 
done by Weinberg in [20]. Goffin, Segre and Weldon [21] have shown that only loops containing 
physical scalars give nonzero contribution to the strong CP phase. Yukawa coupling and mass 
matrix of the up-type quarks are real. So they cannot contribute to strong CP phase. Due to 
this fact the strong CP phase becomes 

5Q = -ImjTr^- 1 !])}. (57) 

where jtft is the 4x4 down-type quark mass matrix as in Eq. ()10|) and X is now the 4x4 
down-type quark self energy matrix. In this model there are three scalar h, s, and t which are 
originated from fluctuation about the vacuum, <pi — ► v + h and S — > -/n{V + s + it) exp(zct). 
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These fields are not mass eigenstates. The physical scalars H k are related to these scalars 
h a (h, s,t) by orthogonal transformations, H k = R ka h a . In this basis the Yukawa couplings T a 
of the Lagrangian change to T k , where T k = RkJ^a- The contribution of the physical scalars to 
the self energy matrix is j22j 



= - y ~7~r~\o I ^[C 1 - x)^V{ + T h jt\ \n{^^x 2 + M?(l - x)\T k . (58 
u 47T Y Jo 



So the strong CP phase will be 

1 n 



5Q 



E 



(4tt) s 



dxlmlTri^Tk^ \n{JZ JZ* x 2 + M 2 (l - x))Y k }\. 



(59) 



We can write the strong CP phase in terms of original Yukawa couplings of the Lagrangian. 
Changing also the logarithmic part of Eq. (J55j) in the diagonal basis the dominant part of the 
strong CP phase comes from D quark mass. Then the strong CP phase will be 



se 



E 

k,a,b 



RkaRkb 

(4vr) 2 



f 1 dxlmV^ ] T a J^- 1 T b J^ ] ^} AA \ii(M D 2 x 2 + M A 2 (1 - x)). (60) 
Jo 



The above expression gives a nonzero value only when a = t and b = h. The explicit calculations 
of strong CP phase have been done earlier jB] . To keep the hierarchy between the scales v and the 
scale V and to avoid the fine tuning in the stationarity equations of the scalar potential in Eq. 
© we should consider two parameters c\ and c 2 of the scalar potential small as 0(v 2 /V 2 ). Light 
scalar mass Mi remains of the order of v. The other two scalar masses become M 2 ~M3~U. So 
the calculation of the order of magnitude of RuRkh^ sin(2a)p^- jHj- In performing integrations 
we keep the assumption M 2j3 >2M^ same as in [6j but 2Mb > Mi is not valid for all f3 as in 
Table 0] and as we consider Mi = Mh = 150 GeV. Under the considerations M 2 and M3 are 
nearly same and M2^^>2Md we get the following expression of strong CP phase 



5Q 



sin 2a 2 _ f/ 2^ 

16?r 2 U ^V 2 



1 + ln 



V 2 



(61) 



where 



h 



Ml 
Ml 

Ml 
Ml 



2 111 { m 1 + I -w 



Mi 



- 1 



AMI 



tan 



AMI 
M'l 


-) 




AM 2 n 
Ml 



for 2M D > Mi, 
for 2M D < Mi. 



(62) 



The integrals are performed by Bento, Branco and Parada in jBj. Using our ansatz of the mass 
matrix of down-type quarks we have strong CP phase of the following form 



5Q 



sin 2a 



in' 



A/(A/ + 2/) 



M 2 D v 2 
V 4 



1 + ln 



' V 2 ' 



(63) 
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p 


Lower Bound 


Upper bound 


Range of 




Ran 


ge of 


Lower Bound 




of Md [GeV] 


of \U bs \ x 10 5 


Am Ba x 10 10 [MeV] 




Ar s ,, 


x 10 2 


of V [TeV] 


— TDD 


oy.i 


n 8^8 


89 tn 1 4fi 4^Q 
oz.uoo lo iiu.^jy 


12.5139 to 23.9984 


1 n 1 

1U. 1 


7 ^ 


79 Q 
( z.y 


n rq^ 


89 D^^ tn 1 4fi 44Q 


12.5145 to 23.9995 


1 9 Q 


— 1{ n 

— o.u 


IUO. 1 


U. JUU 


89 093 tn 1 4fi 4^9 

OZ.UZO LO l^U.^OZ 


12.5155 to 24.0016 


1 7 4 


— ^.0 


1 9^ 4 


u.^±ui 


89 fl91 tn 1 4fi 497 


12.5158 to 24.0022 


1 1 
iy . 1 


4 n 


1 4fi 
l^U.U 


n 41 4 


89 fl1 8 tn 1 4fi 499 
OZ.UlO LO l^U. ^iZZ 


12.5161 to 24.0028 


91 

Zl.U 


— 0.0 


1 87 ^ 


D ^R4 

U.OU^J: 


89 D1 ^ tn 1 4fi 41 7 


12.5165 to 24.0034 


9^ ^ 
zo.o 


— o.z 


9^4 9 
ZO^.Z 


u.ooo 


89 fl1 ^ tn 1 4fi 41 3 
oZ.UlO LO l^U.^ilO 


12.5167 to 24.0038 


9Q n 
zy.u 


9 9 


979 8 
z i z.o 


n ^7Q 
u.o / y 


81 Qfi^ tn 1 4fi ^98 
Ol.yuO LO 14U.OZO 


12.5219 to 24.0140 


ou.o 


9 ^ 




n 491 


81 Qfi^ tn 14fi ^94 
oi.yuo lo i^u. oz^ 


12.5222 to 24.0146 


94 7 

Z^. i 


o.u 


1 ^fi 9 
IOU.Z 


n 4Q^ 
u.^yo 


81 Q^8 tn 1 4fi ^1 fi 
oi.yoo lo iiu.oiu 


12.5226 to 24.0155 


1 Q Q 
ry.y 


0.0 


1 n8 8 

IU0.0 


u. ouo 


81 tn 1 4fi ^flQ 
oi.yoo lo i^u.ouy 


12.5231 to 24.0163 


1 7 1 


4 n 


Q1 Q 

y ± . V 


u.uo J- 


81 Q^1 tn 1 4fi 

Ol. Jjl LU l^iU.OUi 


12.5235 to 24.0171 


1 ^ 1 


4.5 


80.2 


0.697 


81.947 to 146.297 


12.5238 to 24.0178 


13.7 


5.0 


71.6 


0.761 


81.945 to 146.292 


12.5241 to 24.0184 


12.6 


7.5 


49.7 


1.034 


81.937 to 146.278 


12.5249 to 24.0200 


9.3 


10.0 


41.2 


1.220 


81.938 to 146.280 


12.5248 to 24.0197 


6.9 



< 7.204 x 10~ 
Ar Ss 



Table 4: The lower bound on M D using the obtained upper bound of Um, 
for the different values of (3. The upper bound of Uf, s , range of Am Bs , range of 
and the lower bound of the scale V for the different values of j3. We use the solutions 
inputs which are in the rightmost column of Tabled We use aiMz) — 1/128.91, a s [Mz) — 
(0.1187 ±0.002), M z = (91.187 ±0.0021) GeV, M w = (80.425 ± 0.038) GeV, M H = 150 GeV, 



w 



1 - M^/Mf, mt 



(174.3 ± 5.1) GeV and m b = 4.1 to 4.4 GeV p. 



The strong phase will vanish for A/=0, a=0, n/2, 7r, 2tt. Weak phase will also vanish for those 
values of Af and a because for each value and hence become real symmetric. In 

terms of x parameters the form of strong CP phase using the relations in Eq. (|27|) we have 



ah 1 M > 2 




(64) 



Present value of the experimental bound on electric dipole moment of neutron gives value of 
upper bound on strong CP phase, 56 < 2 x 10~ 10 . We can convert the bound of 50 to bound 
of the scale V using the Eq. (p>fj). the solutions of the Eq. (JTHJ) and the bound of in Table 
EJ We show the lower bounds on the scale V for the different values of f3 in Table |3J The 
bounds on the scale V for the different sets of solution for the given inputs will differ within 1 
TeV. There exists a tiny window on the scale V for the bound of 50 for the expression in Eq. 
(J64j) which is small like v. We discard this result because the result come for fine cancellation 
among the terms inside the parenthesis in Eq. (|64p. 
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Upper bound 


Range of 




Ran 


ge of 


aj for 


a n for 




of \U bd \ x 10 5 


Am Bd x 10 10 [MeV] 






x 10 2 


B d J/i/jK s 


B d -> 7T + 7T~ 


1U.U 


3 652 

tj • \JxJ i-i 


2 273 to 5 274 

— . — ( -J u VJ U • I t: 


0.3864 to 0.9108 


7434 


0.2641 


— 7 5 


3 027 

(J.VJZj I 


2 275 to 5 27Q 

— — / ■ J uVJ xj • j-j i U 


0.3863 to 0.9108 


741 8 

VJ. 1 t: X O 


261 5 

VJ . ZjVJ X cj 


— Fi n 

<J . w 


2 204 


2 278 to 5 287 


0.3861 to 0.9102 


7403 

vj. i tiwo 


2597 


—4 Fi 

Tt . cj 


2 00Q 

Zj . VJVJ -J 


2 27Q to 5 28Q 

Zj . Zj 1 ,J L VJ sJ.ZjVJtJ 


0.3859 to 0.9100 


7400 

VJ. 1 tiVJVJ 


2592 

VJ . Zj tj tJ Zj 


—4 n 

TC . VJ 


1 802 


2 280 to 5 2Q1 

Zj.ZUU u VJ sJ.ZjtJX 


0.3859 to 0.9097 


7396 

VJ. 1 tJC/VJ 


2587 

VJ . Zj <JO ( 


— 3 Fi 

*J . tJ 


1 585 


2 282 to 5 2Q4 

Zj.ZjCJZj L VJ sJ.ZjtJt: 


0.3858 to 0.9095 


7392 

VJ. 1 k)U£i 


2581 

VJ . Zj <J(J X 


tj . _ 


1.449 


2 282 to 5 205 

Z >ZUZi UU xj • £j U\J 


0.3857 to 0.9093 


7390 

VJ . 1 '.J ' V 7 


2578 

VJ . Zj tj I (J 


2.2 


1 640 

X . VJ^<J 


2 307 to 5 353 


0.3829 to 0.9026 


7359 

VJ . 1 '.J ' J < J 


2534 

VJ . Z. »J • ) t: 


2 Fi 


1 836 

X . (JOVJ 


2 300 to 5 358 


0.3827 to 0.9022 


7360 

VJ. I (JVJVJ 


2535 


3 n 

o . vj 


2 146 

Zj . X t: VJ 


2 31 2 to 5 365 

Z.iJlZ uU *J.*JVJtJ 


0.3824 to 0.9013 


7361 

VJ . 1 (JVJ X 


2536 

VJ . ZjtJOVJ 


3 Fi 

*J . xJ 


2 451 

Zj ■ t: tJ X 


2 31 6 to 5 373 

Z.U1U uU iJ.O 1 l) 


0.3820 to 0.9005 


7363 

VJ. I (JVJO 


2540 

VJ . ZjfjTiVJ 


4.0 


2.748 


2 310 to 5 381 


0.3817 to 0.8997 


0.7367 


0.2545 


4.5 


3.037 


2.323 to 5.388 


0.3814 to 0.8990 


0.7371 


0.2551 


5.0 


3.316 


2.326 to 5.396 


0.3811 to 0.8982 


0.7377 


0.2559 


7.5 


4.507 


2.338 to 5.423 


0.3798 to 0.8952 


0.7416 


0.2615 


10.0 


5.317 


2.341 to 5.431 


0.3792 to 0.8936 


0.7458 


0.2676 



Table 5: The upper bound of Ubd, range of Amg , range of r Bd x 10 2 , CP asymmetry 

B d 

parameters aj and a ff for the decay channels B® — > J/ipK s and — ¥ 7r+7r_ respectively for the 
different values of (3. We use the solutions for the inputs which are in the rightmost column of 
Table [TJ Mo values are taken from the Table |U to calculate Um bound. We use also the inputs 
in the caption of Table H 0. 



5 Conclusion 



The model of BPR [3: claims that all CP violating phenomena can originate from a single phase 
which appears in the VEV of SM gauge singlet complex scalar S. Bento, Branco and Parada 
started this work where they showed that the weak CP phase and the strong CP phase can 
have a common origin. BPR [3] extend this to the leptonic sector. 

We have made quantitative studies of this model in the quark sector. We find the La- 
grangian parameters for a specific kind of ansatz of down-type quark mass matrix where CKM 
phase is generated minimally. The Lagrangian of the BPR model is CP invariant. CP is bro- 
ken spontaneously through a single phase a in the VEV of the singlet scalar S. Hence, all 
Lagrangian parameters are real and the phases are only due to a. We observe that in the inter- 
val —3.2 < (3 < 2.12 the ansatz of Eq. (J23J) becomes inconsistent. We find the numerical value 
of the Lagrangian parameters outside this interval. We also find that the contribution of the 
bare part to the mass of D quark is negligible compare to the Yukawa part. This observation is 
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independent of f3. One point we should make about the unitarity of Eq. (|12j). To calculate 
U we use the solution of the Eq. (JTSJl where K is unitary. We calculate approximate R, S and 
T of Eq. flU and Eq. (ED) with the solutions Set-1 for = 3.5 in Tableland the value of M D 
for the same solution in Table |U Putting these R, S and T and unitary K in Eq. (fl^jl to test 
its unitarity we see that and ty^ 6 }/ deviate from identity by the amount only 0(1CT 4 ) 

for a few elements. This deviation is not so sensitive to the parameter /3. 

We have observed that this kind of vector quark model cannot increase the partial decay 
width of Z —>■ bb. So we find out the bound on C4& from the condition that the new physics 
decrease the SM value up to experimental lower bound of T Z ->ib considering la error of the total 
decay width of Z and the branching ratio of this decay mode of Z. Here we also should point 
out about value of Mr- Suppose we fix the value of (3 at 3.5 and change Mh from 150 GeV to 
100 GeV. Then the decay width of Z — > bb decreases by nearly 0.01 MeV, remains almost 
same whereas V increases by nearly 0.34 TeV. On the other hand, for the same value of /3, if Mu 
is changed from 150 GeV to 300 GeV, it increases the decay width of Z — > bb by nearly 0.007 
MeV, M B remains almost same whereas V decreases by nearly 0.62 TeV. Scope of new physics 
remains in determining the mass difference Amg, for both q = d,s. Here we have only small 

shift of the Am^ range from SM range. Actually huge uncertainty of \J B Bd f Bd — 221 ± 28±° 2 
and \ / B Bs f Bs = 255 ± 31 introduce huge uncertainty in theoretical value of Am^ and Am Bs 
respectively. This model has very small effect on decay width difference in both B® — B^ and 
B® — B® system. The results are shown in Table El and respectively. We also have observed in 
Table El that the numerical values of the CP asymmetry parameters aj and a n for the decays 
B® — > J/il)K s and B® — ► 7r + 7r~ respectively have effect on the third and somewhere on the 
second place after the decimal from SM value. This part of our analysis is quite general and 
applies to any model containing extra down-type quarks, e.g., models inspired by E 6 grand 
unification. 

The strong CP phase in this model is suppressed by inverse powers of V. It should be noted 
that we have used the solutions of the elements of the quark mass matrix to obtain lower bounds 
on V. There are, however, direct experimental limits on the mass of strongly produced massive 
quarks of charge —1/3. This lower bound of 199 GeV [2H] is derived for a fourth generation 
down-type quark which is produced strongly in pair and decay to bZ via 1-loop FCNC. If this 
bound is assumed to hold for Mp in the BPR model where FCNC exists at the tree level, the 
lower bound on V increases. For example, with (3 = —10, we now get V > 20.3 TeV, whereas 
for /3 = +10, we get V > 15.4 TeV. 

Note added: After completion of our work a new analysis of a group was published [21], 
where the experimental value of m;, has changed to a larger value with a smaller error bar, 
?7i6 — 4.591 ± 0.040 GeV in the kinetic scheme at scale 1 GeV. But we have used the PDG [lj 
value 777.6 — 4.1 to 4.4 in our numerical analysis in this model, and the random value chosen in 
Tabled for finding Lagrangian parameters lies outside the la range advocated in [21]. 
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